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(54) Method for determining the progression of the load torque in an internal-combustion engine 



(57) Method for determining the progression of the 
load torque (Tc(t)) in an internal-combustion engine (1), 
comprising, for each combustion cycle of the engine, 
the steps of detecting (15, 23) the progression of the 
instantaneous angular velocity (v(t)) of the crankshaft 
(6) in at least one given time interval (I) of the combus- 
tion cycle; determining (24,25,26,27) the progression of 
the pressure (P(t)) inside the cylinders (3) of the engine 
(1) for the combustion cycle on the basis of the progres- 
sion of the instantaneous angular velocity (v(t)); and 
determining (19) the progression of the load torque 
(Tc(t)) for the combustion cycle on the basis of the pro- 
gression of the instantaneous angular velocity (v(t)) and 
the pressure (P(t)) inside the cylinders (3). 
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Description 

[0001] The present invention relates to a method for determining the progression of the load torque in an internal- 
combustion engine. 

5 [0002] Methods are known for determining the progression of the load torque applied to an internal-combustion 
engine during a combustion cycle of the engine, whereby the progression of the load torque is obtained both as a func- 
tion of the progression of the angular velocity of the engine during the combustion cycle itself and as a function of the 
progression of the pressure inside at least one cylinder. 

[0003] In particular, these methods envisage the use of at least one sensor outside the engine block so as to be able 
10 to reconstruct the progression of the pressure inside the cylinder. In fact, in some cases, a pressure sensor is mounted 
directly inside the combustion chamber of the cylinder so as to output a signal indicating the internal pressure. In other 
cases, the pressure inside the cylinder is obtained by means of a piezoelectric accelerometer mounted on the engine 
block, opposite the cylinder, so as to output an electric signal representing the vibrations of the engine block, which, as 
is known, are directly correlated to the progression of the internal pressure. 
15 [0004] Therefore, the said methods for determining the load torque require a sensor outside the engine block, which 
generally involves not insignificant costs due to the sensor as such, the installation operations and any modifications 
which must be performed to the engine block in order to provide the seats suitable for receiving said sensor. Moreover, 
both the pressure sensor and the accelerometer, in view of the hostile conditions of the environment in which they oper- 
ate, are subject to continuous mechanical and/or thermal stresses which could adversely affect the correct functioning 
20 thereof as well as reduce their average working life. 

[0005] The object of the present invention is that of providing a method for determining the progression of the load 
torque in an internal-combustion engine, which does not have the drawback described above and by means of which it 
is possible to determine the load torque without modifying the engine block and without the aid of sensors outside the 
engine block itself. 

25 [0006] According to the present invention, a method for determining the load torque in an internal-combustion engine 
is provided, characterized in that it comprises the steps of detecting the progression of a physical parameter indicating 
the instantaneous angular velocity of the said engine; determining the progression of the pressure inside at least one 
cylinder of the said engine as a function of the progression of the said physical parameter; and determining the progres- 
sion of the said load torque as a function of the progression of the said physical parameter and the progression of the 

30 said internal pressure. 

[0007] The present invention will now be described with reference to the accompanying drawings, which illustrate a 
non-limiting example of embodiment thereof and in which: 

Figure 1 shows in schematic form, with parts removed for the sake of clarity, an internal-combustion engine pro- 
35 vided with a processing device which implements the method according to the present invention; 

Figure 2 shows in schematic form a block diagram of the operations performed by a first circuit forming part of the 
processing device according to Figure 1 ; 

Figure 3 shows a cartesian diagram representing the operating points of the engine according to Figure 1 ; 
Figure 4 shows the progression of the pressure inside a cylinder of the engine as a function of the angular position 
40 of the crankshaft during a combustion cycle; 

Figure 5 shows in schematic form a second circuit forming part of the processing device according to Figure 1 ; 
Figure 6 shows the progression of the torque transmitted to the crankshaft as a function of the angular position of 
the crankshaft during a combustion cycle; 

Figure 7 shows a circuit diagram of a block forming part of the second circuit according to Figure 5; and 
45 - Figure 8 shows a progression of the load torque reconstructed in accordance with the method according to the 
present invention. 

[0008] With reference to Figure 1, 1 denotes in its entirety a "four-stroke" internal-combustion engine of the known 
type, comprising a support base 2 and a plurality of cylinders 3 (only one of which is shown) which are supported by 
so the base 2 and define, at their top ends, respective combustion chambers 4. Each cylinder 3 has, mounted internally, a 
respective piston 5 which is connected to the crankshaft 6 by means of a respective connecting rod 7 in a known man- 
ner and is movable axially inside the cylinder 3 itself between an upper end-of-travel position (known as "top dead cen- 
tre") and lower end-of-travel position (known as "bottom dead centre"). 

[0009] The engine 1 comprises an intake manifold 8 connected to the cylinders 3 in order to supply a flow of the com- 
55 bustion agent (air) into the cylinders 3 themselves, a device 9 for supplying the fuel into the cylinders 3, an ignition 
device 12 for triggering combustion inside the chambers 4 and an exhaust manifold 13 connected to the cylinders 3 for 
conveying away the combusted gases. 

[0010] The engine 1 is controlled by a central control unit 14 co-operating with a plurality of sensors inside the engine 
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block so as to receive at its input data signals representing physical parameters of the engine 1 such as, for example, 
the position of the butterfly valve, the temperature of the cooling liquid, and so on. In particular, in accordance with that 
shown in Figure 1 , the central control unit 1 4 is connected to a sensor 1 5 for the velocity of the crankshaft 6 which out- 
puts a velocity signal v(t) representing the instantaneous angular velocity of the shaft 6 and to a pressure sensor 16 
s arranged along the intake duct 8 so as to produce a pressure signal Pc(t) representing the pressure inside the manifold 
8 itself. 

[001 1 ] According to the present invention, the central control unit 14 comprises a processing device 1 7 which receives 
at its input the instantaneous angular velocity signal v(t) of the shaft 6 and is able to output a signal Tc(t) representing 
the load torque applied to the engine for each combustion cycle of the engine. 

w [001 2] The processing device 1 7 comprises a circuit 1 8 which receives at its input the angular velocity signal v(t) and 
is able to output a signal P(t) indicating the pressure inside the cylinder 3. The device 1 7 comprises, moreover, a circuit 
19 which receives at its input the signal v(t) and the signal P(t) and is able to output the signal Tc(t). 
[0013] Below, with reference to Figure 2, the circuit 18 which is able to reconstruct the progression of the pressure 
P(t) inside the cylinder 3 during a combustion cycle of the engine will be described. 

is [0014] With reference to Figure 2, an initial START block is followed by a block 22 in which the existing operating point 
of the engine 1 is determined in relation to a set of operating parameters of the engine 1 itself. 
[0015] In the example illustrated the existing operating point is determined by two operating parameters: the rpm of 
the engine (number of revolutions per minute) and the mean pressure Pern inside the intake manifold 8. In particular, 
the existing operating point is determined on the basis of the values of the mean pressure Pern and the rpm of the 

20 engine relating to the combustion cycle prior to the existing cycle, said values being available in that they are calculated 
by averaging the signal Pc(t) and the signal v(t). 

[001 6] It is obvious that it is possible to determine the existing operating point in relation to a different number of oper- 
ating parameters or simply in relation to a different pair of parameters themselves. 

[0017] The block 22 is followed by a block 23 in which a given number N of values of the signal v(t) indicating the 

25 instantaneous angular velocity are sampled; these sampled values, which are indicated below by v 0 y^ t V N . 1( relate 

to a time interval I which is inside the combustion cycle and which comprises the instant at which the piston 5 is located 
at the top dead centre during ignition. In the example shown, the interval I corresponds to that time interval in which the 
angular position 6 of the shaft 6 (Figure 1) lies within an angular interval l 9 (Figure 4) having an amplitude equal to 180° 
and centred in the angular position in which the piston 5 is located at the top dead centre. 
30 [0018] The block 23 is followed by the block 24 in which the discrete Fourier transform v(jnm) of the angular velocity 
signal v(t) "sampled" (or considered) in the interval I is calculated, in accordance with the expression: 



N-1 

I 

35 n-0 



V(A>J=£ v n .e" /( " mn) (1) 



where nm is the angular frequency defined as 

n m = fj— ( 2 ) 

where N is the abovementioned number of sample values of the signal v(t) sampled in the interval I and m is the 

45 number of the m-th harmonic present in the signal v(t) itself. 

[0019] The block 24 is followed by a block 25 in which, as will be explained below, a frequential response function 
H(jftm) relating to the existing operating point of the engine I and expressing the relationship between the progression 
of the angular velocity signal v(t) and the pressure signal p(t) inside the cylinder 3 in the domain of the angular frequen- 
cies is calculated. As will be explained more clearly below, a given frequential response function H(jOm) is associated 

so with each operating point of the engine. 

[0020] The block 25 is followed by the block 26 in which, based on the knowledge of the Fourier transform V(jOm) and 
the frequential response function H(jnm), the discrete Fourier transform Pflflm) of the signal P(t) representing the pres- 
sure inside the cylinder 3 is calculated, in accordance with the expression: 



55 



4 



8/18/2007, EAST Version: 2.1.0.14 



EP0 985 919 A1 



[0021 ] The block 26 is followed by a block 27 in which the anti-transformation operation in accordance with the Fourier 
transform P(jnm) is performed so as to obtain the evolution of the signal P(t) for the pressure inside the cylinder 3 in the 
time interval I. In particular the anti-transformation operation is performed in accordance with the known expression: 

Pn-JZ P{&m)-e mm ' n) **hn=0 AM (4) 

where P 0 .Pi,- .. .Pn-1 indicate the reconstructed values of the pressure inside the cylinder 3 at the instants when 

10 the values v^V^ V N .-, within the time interval I were sampled. 

[0022] In this way, by means of the signal P(t), and in particular by means of the values P 0 ,Pi Pn-l the circuit 18 

(Figure 1 ) provides at its output the progression of the pressure P(t) inside the cylinder 3 for the time interval I and, basi- 
cally, for the entire combustion cycle, the progression of the internal pressure itself outside the interval I during the 
intake and exhaust stages being known a priori. 

15 [0023] In accordance with that described above, the relation (3) is therefore the key for being able to reconstruct the 
internal pressure signal P(t) from the angular velocity signal v(t) output by the sensor 15. The relation (3) is based on 
the existence of a linear correlation, which has been proven experimentally, between the progression of the instantane- 
ous angular velocity v(t) and the pressure p(t) inside a cylinder 3 in the interval I. Experimentally it has been noted that, 
in order to reconstruct the signal P(t) with a reasonable degree of approximation, rt is possible to take into consideration 

20 only a certain number of harmonics, for example the first seven harmonics; as a result the expression (4), by means of 
which it is possible to determine the progression of the internal pressure, becomes: 

6 

Pn"J I P(jnj-e nnm ' n) withn=0 AM. 

25 m=0 



[0024] The calculation of the function H(jOm), which is carried out in the block 25, is performed on the basis of the 
knowledge of a plurality of frequential response reference functions H R (jOm), each of which relates to a respective ref- 

30 erence operating point R of the engine 1 and is stored inside the processing device 17. 

[0025] In particular, by representing the operating points R in the plane (rpm, Pern) of the possible operating points 
of the engine (Figure 3), the function H(jOm) relating to the existing operating point (indicated by K in the plane (rpm, 
Pern)) is calculated by taking into consideration a subgroup S of the points R. This subgroup S comprises the reference 
points R closest to the existing operating point K. 

35 [0026] The function HfjOm) relating to the point K is calculated by determining a weighted mean of the reference func- 
tions H R Qnm) relating to the points R belonging to the subgroup S, in accordance with the expression 



S H R (jCl m )-D R (5) 

40 ReS 

where D R is the weight associated with the point R and is a function of the distance between the point R itself 
and the point K in the plane (rpm, Pern). 
45 [0027] In particular, in the example illustrated in Figure 3, the reference points R are arranged so as to form a grid Q 
and the subgroup S comprises the four points R which surround the existing operating point K. 
[0028] To summarise, the calculation of the function H(jfim) relating to the existing operating point of the engine is 
thus performed (Figure 3) as follows: 

so - the existing operating point K is determined in the plane (rpm, Pern) (block 22) ; 

the subgroup S is determined by defining the reference points R closest to the point K, 

- the weights D R are calculated as a function of euclidean distances between these points R and the point K in the 
plane (rpm-Pcm); and 

- the function Hfjfim) is calculated on the basis of the reference functions H R (jOm) associated with the points R of 
55 the subgroup S in accordance with the expression (5). 

[0029] From the above it can therefore be seen that, on the basis of the knowledge of the frequential response func- 
tions H R (jQm) at the points R in the grid G, it is possible to determine the progression of the pressure P(t) inside the 
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cylinder 3 during the combustion cycle, once the progression of the angular velocity v(t) in the interval I is known. 
[0030] The reference functions H R (jOm) at the reference points R, stored in the processing device 1 7, each express 
the relationship between the Fourier transform VQnm) of the angular velocity signal v(t) and the Fourier transform 
Prjfim) of the internal pressure signal P(t) when the engine 1 is at the operating point R, namely: 

5 

w [0031] These Sequential response functions H R (jOm) were obtained with reference to a sample engine having the 
same geometric and structural characteristics as the engine 1 (for example connecting-rod length, crank length, piston 
area, etc.). In fact, such a sample engine was equipped with a pressure sensor mounted inside the combustion cham- 
ber of the cylinder so as to produce directly at its output the signal P(t) for the pressure inside the cylinder. Once the 
operating parameters of the sample engine have been fixed so that the associated operating point coincides with the 

is reference operating point R, the progressions of the instantaneous angular velocity v(t) and the internal pressure P(t) 
for the observation interval I are detected directly in the sample engine. At this point the Fourier transforms V(jftm) and 
P(jOm) of the angular velocity v(t) and the internal pressure P(t), respectively, are calculated and then the function 
H R (jflm) relating to the point R is obtained simply from the relation: 



20 



[0032] Then, all the Sequential response reference functions H R (jnm) are determined with reference to the sample 
25 engine. These functions H R (jOm) are stored in the processing device 17 of the engine 1 so that it is possible to deter- 
mine the temporal progression of the pressure inside the cylinder 3 without installing any pressure sensor inside the 
combustion chamber 4. 

[0033] With reference to Figure 5 the circuit 19 (Figure 1) able to determine the progression of the load torque Tc(t) 
from the angular velocity v(t) (measured by means of the sensor 15) and from the internal pressure P(t) (determined by 
30 means of the circuit 18) will now be described. 

[0034] The internal pressure signal P(t) is supplied to.a conversion circuit 30 of the known type, which is able to output 
a signal Tt(t) indicating the torque transmitted to the crankshaft 6. This conversion circuit 30 is able to determine the 
progression of the transmitted torque Tt(t) by applying the known law of the crank mechanism to the internal pressure 
P(t). 

35 [0035] The angular velocity signal v(t), on the other hand, is supplied to an adder input of an adder node 31 which 
has, moreover, a subtracter input connected to the output 32u of a block 32 which models in a known manner the 
dynamic behaviour of the engine-load system from the rotational point of view. In particular, the block 32 has two inputs 
32a and 32b, the input 32a of which receives the transmitted torque signal Tt(t) output by the conversion circuit 30, while 
the input 32b receives a signal Tc 8t (t) which, as will emerge clearly below, represents the estimated torque load. The 

40 block 32 is able to generate at its output 32u a signal V 8t (t) indicating the estimated angular velocity for the shaft 6 on 
the basis of the values assumed by the signals Tt(t) for the transmitted torque and Tc 8t (t) for the estimated load torque. 
[0036] The node 31 is able to provide at its output a signal Av(t) defined by the difference between the signal v(t) out- 
put by the sensor 15 and the signal v 8t (t) indicating the angular velocity estimated on the basis of the engine-load 
model. This signal Av(t) represents the error in angular velocity which exists between the angular velocity v actually 

45 measured and the angular velocity v 8t estimated on the basis of the reconstructed transmitted torque Tt(t) and the esti- 
mated load torque Tc 8t (t). 

[0037] The error signal Av is supplied to an estimator circuit 33 which is able to process it so as to provide the esti- 
mated load torque signal TCg, at its output 33u. The output 33u is fed back to the input 32b of the block 32 so as to allow 
the block 32 itself to receive the estimated load torque signal Tc 8t . 
so [0038] The estimator circuit 33 (described in detail below) is designed to operate in such a way that the estimate of 
the load torque Tc^ is changed so that the estimated angular velocity V 8t tends to reach the measured angular velocity 
v. 

[0039] The output 33u, moreover, is connected to a low-pass filter 34 designed to eliminate high-frequency spectral 
components from the signal Tc^ so as to output the signal Tc, namely the reconstructed load torque which is actually 
55 applied to the engine. 

[0040] Figure 7 shows the circuit diagram of the estimator circuit 33 which is formed using the sliding mode technique 
and, in the example illustrated, is of the proportional-integral type. In particular, the circuit 33 has a multiplier block 36 
which receives at its input the signal Av, is able to multiply said signal by a parameter C 1 which can be calibrated and 
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is connected at its output to an adder node 37. 

[0041] Moreover, the circuit 33 has two blocks 38 and 39 which are arranged in series and the block 38 of which 
receives at its input the signal Av and provides at its output the "sign" of the difference between the measured velocity 
signal v and the estimated angular velocity signal v 8t . The block 39, on the other hand, is a multiplier block which is able 
s to multiply the output of the block 38 by a parameter C 2 which can be calibrated, before supplying it to the adder node 
37. 

[0042] Finally, the circuit 33 has an integrator 40 and a multiplier 41 which are arranged in series and the integrator 
40 of which receives at its input the output of the block 38, while the multiplier 41 is able to multiply the output signal of 
the integrator 40 itself by a parameter C 3 which can be calibrated and then supply it to the adder node 37. The output 
10 of the node 37 defines the output 33u of the estimator circuit 33 and provides the estimated load torque Tc 8t which is 
then calculated in accordance with the expression: 

Tc 8t = C y • Av+ C 2 • s/gn(Av)+ C 3 • Isign(Av) 

15 where "sign" is the sign function belonging to the block 38. 

[0043] Below, operation of the circuit 1 9 which is able to reconstruct the load torque Tc(t) from the measured angular 
velocity v(t) and the internal pressure P(t) will be described. 

[0044] Once the reconstructed internal pressure P is known, the transmitted torque Tt is determined by means of the 
conversion circuit 30 and it is supplied to the block 32 which models the engine-load system. 
20 [0045] At this point the input 32b of the block 32 provides a first trial load torque Tc which is, for example, calculated 
as an offset value with respect to the reconstructed transmitted torque Tt. The block 32, in response to these inputs, 
provides at its output a first estimate of the angular velocity, namely the angular velocity which the shaft 6 should have 
in order to balance the reconstructed transmitted torque Tc with the trial load torque Tc. 

[0046] The adder node 31 then calculates the velocity error Av existing between the measured velocity and the esti- 
25 mated velocity and supplies this error to the estimator circuit 33. The circuit 33 in turn processes this velocity error Av 
so as to output a second trial load torque which will be fed back to the input 32b of the block 32. This feedback will result 
in correction of the estimate of the angular velocity output from the block 32. At this point, the operations described 
above are repeated until the velocity error Av becomes approximately zero. 

[0047] When the estimated angular velocity (output from the block 32) is equal to the measured angular velocity, the 
30 estimated load torque output from the circuit 33 is approximately equal to the load torque which is actually applied to 
the engine and, consequently, the process for determination of the load torque ends. 

[0048] It should be emphasized that the rapid convergence of the process which results in the load torque being 
obtained (namely the rapid convergence, to zero, of the velocity error Av) is ensured by the estimator circuit 33 (Figure 
7), the control action of which, as already specified above, is expressed by the relation: 

35 

Tc at = C, • Av+ C 2 • sign{Av)+ C 3 • I sign(Av) 

[0049] In fact, the term • Av has the function of forcing the load torque estimate to change so as to nullify the veloc- 
ity error Av as quickly as possible, while the terms C 2 • sign(Av) and C 3 • J sign (Av) associated with the sign of the 

40 error Av perform a control action aimed at ensuring that, after a certain number of oscillations about the final value of 
the estimate for the load torque, the estimated load torque coincides with the final value itself. 
[0050] From what has been described above it can be seen that in the engine 1 , which is intended to be a mass-pro- 
duced vehicle engine, it is possible to determine the progression of the load torque applied to the engine, cycle after 
cycle, on the basis solely of the measurement of the angular velocity of the crankshaft and without the aid of any sensor 

45 on the outside of the engine block such as, for example, a pressure sensor mounted directly inside the combustion 
chamber or a piezoelectric accelerometer. 

[0051 ] In accordance with that illustrated with reference to Figures 4, 6 and 8, by means of the method according to 
the present invention it is possible to reconstruct the progressions of the pressure inside the cylinder (Figure 4), the 
torque transmitted to the crankshaft 6 (Figure 6) and the load torque (Figure 8) with an excellent degree of approxima- 

50 tion. In fact, representing graphically as a continuous line the progressions of the internal pressure P and the transmit- 
ted torque Tt as a function of the angular position e of the shaft 6, it is possible to note that these progressions do not 
deviate significantly from the progressions which are obtained by measuring directly these parameters (shown in bro- 
ken lines in Figures 4 and 6). Similarly, with reference to Figure 8, representing graphically as a continuous line the 
reconstructed progression of the load torque Tc as a function of time, it is possible to note that this progression does 

55 not deviate significantly from the actual progression which is obtained by measuring the load torque (shown in broken 
lines). 

[0052] As a result of the method described it is therefore possible to determine the evolution of the internal pressure, 
the transmitted torque and the load torque without having to incur costs due to the installation of sensors outside the 
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engine block itself and on the basis solely of measurements of the angular velocity of the crankshaft. This ensures that 
the method can be applied to mass-produced engines in a simple and low-cost manner. 

Claims 

5 

1. Method for determining the progression of the load torque in an internal-combustion engine (1), characterized in 
that it comprises the steps of: 

a) detecting the progression of a physical parameter indicating the instantaneous angular velocity (v(t)) of the 
10 said engine (1); 

b) determining the progression of the pressure (P(t)) inside at least one cylinder of the said engine as a func- 
tion of the progression of the said physical parameter; and 

c) determining the progression of the said load torque (Tc(t)) as a function of the progression of the said phys- 
ical parameter and the progression of the said internal pressure (P(t)). 

15 

2. Method according to Claim 1 , characterized in that the said physical parameter is the instantaneous angular veloc- 
ity (v(t)) of the crankshaft (6). 

3. Method according to Claim 2, characterized in that the said step of determining the internal pressure (P(t)) accord- 
20 ing to paragraph b) comprises, for each combustion cycle of the engine (1), the sub-steps of: 

c) determining (22) the existing operating point (K) of the engine (1) as a function of a set of operating param- 
eters (rpm, Pern) of the engine (1); 

d) detecting (23) the progression of the instantaneous angular velocity (v(t)) of the crankshaft (6); 

25 e) determining (25), with regard to the existing operating point (K), a transfer function (H(jnm)) for the said 

pressure (P(t)) inside the cylinder (3) and the said instantaneous angular velocity (v(t)); and 

f) obtaining (24,26,27) the progression of the pressure (P(t)) inside the said cylinder (3) from the progression 
of the instantaneous angular velocity (v(t)) and from the transfer function (HfJOm)) relating to the existing oper- 
ating point (K). 

30 

4. Method according to Claim 3, characterized in that the said transfer function (H(jnm)) is a frequential response 
function for the said internal pressure (P(t)) and instantaneous angular velocity (v(t)). 

5. Method according to Claims 3 or 4, characterized in that the said step of detecting (23) the progression of the 
35 instantaneous angular velocity (v(t)) according to paragraph d) is performed by sampling (23), for each combustion 

cycle of the engine (1), the said instantaneous angular velocity (v(t)) at a plurality of sampling instants, obtaining a 
plurality of measured values (vq^ V N ..|). 

6. Method according to Claim 5, characterized in that the said substep of obtaining (24,26,27) the progression of the 
40 internal pressure (P(t)) according to paragraph f) comprises the further substeps of: 

g) calculating (24) the discrete Fourier transform of the angular velocity (VQftm)) on the basis of the plurality of 
measured values (v 0 ,v 1 V N .i ; 

h) applying (26) the transfer function to the said Fourier transform of the angular velocity (V(jOm)) so as to 
45 obtain the Fourier transform of the pressure (PGftm)) inside the cylinder (3); 

i) calculating (27) the Fourier anti-transform of the Fourier transform of the internal pressure (PQftm)) so as to 

obtain a plurality of estimated values of the internal pressure (P 0 ,Pi Pn-i) at the said instants for sampling 

of the instantaneous angular velocity (v(t)), the said plurality of estimated values of the internal pressure (P 0 , 
P 1 P N .-,) defining the progression of the pressure (P(t)) inside the cylinder (3). 

so 

7. Method according to Claim 5 or Claim 6, characterized in that sampling of the said instantaneous angular velocity 
(v(t)) is performed, for each combustion cycle of the engine (1), in a time interval (I) during which the piston (5) 

• associated with the cylinder (3) moves in a predefined vicinity of the top dead centre during the ignition phase of 
the engine (1). 

55 

8. Method according to any one of Claims 3 to 7, characterized by the fact of comprising the further preliminary steps 
of: 
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I) determining a plurality of reference operating points (R) of the engine (1), each of which is a function of the 
said operating parameters (rpm, Pern) of the engine (1); and 

m) determining for each reference operating point (R) a respective reference transfer function (H R (jOm)) for the 
pressure (P(t)) inside the cylinder (3) and the instantaneous angular velocity (v(t)) of the crankshaft (6). 

5 

9. Method according to Claim 8, characterized in that the said operating parameters (rpm, Pern) of the engine (1) 
comprise the mean angular velocity (rpm) of the crankshaft (6) and the mean pressure (Pern) in the intake manifold 
(8) of the engine (1). 

w 10. Method according to Claim 9, characterized in that the mean angular velocity (rpm) of the crankshaft (6) and the 
mean pressure (Pern) in the intake manifold (8) define a plane (rpm-Pcm) of possible operating points of the engine 
(1), each said operating point of the engine (1) being represented by a respective point in the plane (rpm-Pcm), the 
reference operating points (R) of the engine (1) defining a grid (G) in the said plane (rpm-Pcm). 

is 11. Method according to any one of Claims 8 to 10, characterized in that the said transfer function (HQOm)) for the 
internal pressure (P(t)) and the said instantaneous angular velocity (v(t)) of the crankshaft (6) at the existing oper- 
ating point (K) of the engine (1) is calculated by processing the said reference transfer functions (H R (jnm)). 

1 2. Method according to Claims 1 0 and 1 1 , characterized by the fact of defining on the said grid (G) a subgroup (S) of 
20 the reference operating points (R) comprising the reference operating points (R) which are closest to the existing 

operating point (K) in the said plane (rpm-Pcm), the transfer function (H(jOm)) at the existing operating point (K) 
being calculated as a weighted mean of the reference transfer functions (H R (jnm)) relating to the reference oper- 
ating points (R) of the said subgroup (S), each reference transfer function (H R (jftm)) acting in the weighted mean 
with a weight which is a function of the euclidean distance between the corresponding reference operating point (R) 
25 and the existing operating point (K). 

13. Method according to any one of the preceding claims, characterized in that the said step of determining the pro- 
gression of the said load torque (Tc(t)) as a function of the progression of the instantaneous angular velocity (v(t)) 
and the internal pressure (P(t)) comprises the substep of obtaining (30) the progression of the torque (Tt(t)) trans- 

30 mitted to the crankshaft (6) from the progression of the pressure (P(t)) inside the cylinder (3). 

14. Method according to Claim 13, characterized in that it comprises, moreover, the steps of: 

n) assigning a trial load torque (Tc 8t ) to a first input (32b) of a block (32) which models the dynamic behaviour 
35 of the system defined by the combination of engine and load, said block (32) having a second input (32a) 

receiving the said torque (Tc) transmitted to the crankshaft; 

o) estimating (32), by means of the said block (32), the angular velocity of the crankshaft (V 8t ) as a function of 
the trial load torque (Tc^ and the torque (Tt) transmitted to the crankshaft (6). the estimated angular velocity 
being the velocity at which the crankshaft (6) should travel in order to balance the transmitted torque (Tt) with 
40 the trial load torque (TCgt); 

p) comparing (31) the measured angular velocity (v) with the estimated angular velocity (v 8t ), generating a 
velocity error signal (Av); 

q) correcting (33) the trial load torque (Tc 8t ) by means of an estimator circuit (33) which receives at its input the 
said velocity error signal (Av) and is able to estimate the load torque {T^ on the basis of the velocity error sig- 
45 nal (Av) itself; and 

r) repeating the steps according to paragraphs o), p) and q) until the velocity error signal (Av) is approximately 
equal to zero so that the trial load torque (Tc^O output by the estimator circuit (33) is approximately equal to the 
load torque actually applied to the engine. 

so 1 5. Method according to Claim 1 4, characterized in that it comprises, in addition, the step of filtering (34) the estimated 
load torque (Tc 8t ) output by the said estimator circuit (33) in order to eliminate high-frequency spectral components. 

16. Method according to Claim 14 or Claim 15, characterized in that the said step according to paragraph q) is per- 
formed by processing the said velocity error signal (Av) proportionally (36,39) and integrally (41) by means of the 

55 said estimator circuit (33). 

17. Method according to Claim 16, characterized in that the estimated load torque (Tc 8t ) is calculated by the estimator 
circuit (33) as a function of the velocity error signal (Av) in accordance with the expression: 
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Tc st = C 1 • Av + C 2 • sign(Av)+ C 3 • I sign(Av) 

where sign(Av) represents the sign function of the velocity error AV and , C 2 and C 3 are parameters which 
can be calibrated. 
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Calculation of the DFT of V 0 ...V N ., 
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Calculation of the frequential response 
function HQClm) at the existing operating 
point 



Calculation of the DFT of the pressure 
inside the cylinder 

P(jn^=Xl£mL 
H,(jn m ) 



Calculation of the IDFT of PQCim) 
P(jA m ) IDfT. P 0 ^....P^ 
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